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ABSTRACT 


The  propagation  of  a  simple  acoustic  transient  in  an  isovelocity 
water  layer  was  investigated  for  the.  purpose  of  studying  the  region 
of  validity  of  a  suggested  correlation  between  exact  solutions  in 
ducts  and  solutions  for  propagation  in  a  layer.   The  Laplace  trans- 
form method  was  used  to  obtain  approximate  solutions  for  the  acoustic 
pressure,  particle  velocity,  and  particle  displacement  resulting 
from  a  step- function  input  in  velocity.   The  radiation  impedance 
characteristics  of  a  Mylar  transducer  were  investigated  to  study 
their  effect  on  the  received  waveforms.   Computer  programs  were 
written  to  evaluate  and  graph  the  resulting  waveforms.   The  wave- 
forms resulting  from  a  step- function  input  in  velocity  were  observed 
and  compared  with  the  predicted  waveforms.   With  the  use  of  a  family 
of  Mylar  transducers,  a  somewhat  doubtful  correlation  was  obtained 
in  the  region  of  validity  of  the  theoretical  solutions.   Possible 
causes  of  the  disagreement  between  the  theory  and  the  experiment 
were  investigated.   It  was  found  that  the  theory  is  believed  to  be 
correct,  that  the  experiment  was  not  accurate  enough,  and  that 
further  investigations  should  be  conducted  to  study  the  assumed 
pressure-release  characteristics  of  \   inch  polyethylene. 
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1.   INTRODUCTION 

Although  most  of  the  current  investigations  of  the  propagation 

•  a  •       (H-13) 

of  energy  in  waveguides  have  been  done  for  electromagnetic  cases 

some  have  considered  the  propagation  of  pulsed  acoustic  energy  in 
waveguides     .   The  dispersion  of  an  ultrasonic  pulse  propagated 
through  a  waveguide  of  circular  cross-section  filled  with  water 
has  been  used  for  the  approximate  determination  of  the  frequency 
spectrum  of  the  original  pulse    .   J.  M.  Proud  et  al.    have  shown 
that  the  dispersion  of  a  rectangular  pulse  along  an  acoustic  wave- 
guide of  rectangular  cross-section  can  be  used  to  calibrate  the  wave- 
guide, so  that  it  can  be  used  to  predict  the  effects  of  similar- 
dispersion  on  the  transmission  of  pulses  of  other  shapes  that  are 
not  amenable  to  actual  computation.   It  has  been  shown,  by  considering 
pulse  compression  in  an  acoustic  waveguide,  that  a  frequency  modulated 

input  pulse  can  be  compressed  into  an  output  pulse  of  shorter  duration 

(3) 
and  of  higher  peak  amplitude    .   In  addition,  studies  of  transient 

(  (\    7  "^ 

propagation  in  layered  media  have  been  made   '   .   Many  of  these 
studies  have  used  approximation  techniques,  such  as  the  method  of 

stationary  phase,  to  solve  the  transient  propagation  problem.   J.  R. 

(8) 
Wait    has  presented  a  summary  of  these  techniques  useful  for  cal- 
culating the  resultant  propagated  waveform.   Using  the  Laplace  trans- 

(5  9  10) 
form  technique,  others     '    have  investigated  transient  radiation 

and  scattering  from  cylindrical  bodies. 

Recently,  exact  solutions  for  the  propagation  of  simple  transients 

in  waveguides  have  been  published   '       .It  has  been  suggested 

that  there  is  a  correlation  within  a  certain  region  of  validity  between 


13 


exact  solutions  in  ducts  and  solutions  for  the  propagation  of 

.   .     ,   .    ,      (5,14) 
acoustic  transients  in  isovelocity  layers 

The  purpose  of  this  research  is  to  extend  investigation  of 

the  region  of  validity  of  the  suggested  correlation  and  of  the 

restrictions  under  which  the  solutions  describing  the  transient 

response  of  the  layer  are  valid. 
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2.   THEORETICAL  DEVELOPMENT 

Development  of  Solutions 
The  response  of  an  isovelocity  layer  to  simple  transient 
excitations  will  be  developed  using  the  following  assumptions: 

(a)  cylindrical  velocity-source  of  radius  a, 

(b)  constant  layer  depth  d,  and 

(c)  pressure-release  upper  and  lower  boundaries. 

The  wave  equation  for  acoustic  velocity-potential  in  cylindrical 
coordinates,  under  the  assumption  of  radially- symmetric  motion,  is 


2.1 


where  r  is  the  radial  direction  and  z  is  the  transverse  direction. 


The  convenient  separation  of  variables 


2.2 


casts  Eq.  2.1  into 

where  the  radial  component  of  the  particle  velocity  becomes 
and  the  acoustic  pressure  becomes 

For  pressure-release  top  and  bottom,  we  have  p  =  0  at  z  =  0 

WnZ  nTTc 

and  z  =  d;  therefore,  Z(z)  =  D  sin  where  w  =  — —   .   Thus, 

c  n     d 

the  equation  for  Z(z),  which  gives  the  dependence  of  the  acoustic 
variables  on  the  layer  depth,  determines  the  eigenfrequencies  w  . 
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These  are  the  cut-off  frequencies  of  the  various  radially-symmetric 
normal  modes  of  the  layer. 

Use  of  the  Laplace  transform  method  to  solve  Eq.  2.2  with 
the  boundary  condition  at  the  source  written  as 

urU,t|=  lU)ML^j 

which    implies   excitation  of   only   one  mode,    casts  Eq.    2.2   into 
the   simpler   differential  equation 


f     ll         ii       / 10*1  +A*-\1  -2  /       v      rs 


2.3 


where^(r,s)  is  the  Laplace  transform  of0(r,t).   V(t)  Z(z)  is 

the  radial  velocity  at  the  source.   The  relevant  boundary  conditions 

become 

¥\     --VM 

and 

$  -*>0   o-s   r  -*•©©  . 

r     ri  2> 

The  change  of  variable  u  =-  J  w  +  s   casts  Eq.  2.3  into  the 


form 


aul  '   U   c*"     -  2.4 


which  has  solution 


<Z>--  AloM+  BKoM 


where  I  and  K  are  the  modified  Bessel  functions.   Since  I  (u)^OO 
00  o 

and  K  (u)— ^0  as  u-vo©,  the  solution  to  Eq.  2.4  which  obeys  the 


required  boundary  conditions  is 


where 


C(6.V-V(0-j===r   MlfTw7)  . 
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Under  the  Laplace  transform,  the  acoustic  pressure,  radial  component 
of  particle  velocity,  and  radial  component  of  particle  displacement 
become 


-V(^ 


Ko(£J^+um£) 


) 


Formal  inversion  results  in  expressions  like 
±ioo 


Ur  = 


-too 

where  the  contour  is  indented  into  the  right  half-plane  around  the 
singularities  on  the  imaginary  axis.   Exact  solutions  to  Eq.  2.5 
do  not  appear  to  be  known,  but  certain  approximations  may  be  used 
to  obtain  approximate  solutions  for  the  leading  portions  of  the 
propagating  transient. 

For  large  values  of  argument,  use  of  standard  approximations 
for  the  modified  Bessel  functions  results  in  the  following  asymptotic 
form 

-  la   e        ;  i  +•  Aw*  (£5  J 

*•  2.6 


(15) 


+  B>/v^(flT5")  *  Cvn*A  ("STff)  *"*A 
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where    <S  -J**M  +  &*"       and  T    1     ILlfl 


and 


Boi=    r 


^  8   *     8     n.  J 


11         lift     "  TTfc  a.       a<3  ^  /w    • 


The  symbol  T  represents  the  time  of  flight  over  the  shortest 
path  between  the  surface  of  the  source  and  the  field  point.   The 
restriction  on  s  necessary  to  satisfy  the  approximation  which  led 
to  Eq.  2.6  is 

Imaginary  s  greater  than  w  in  magnitude  requires 

|*(:*W 

where  , — - — — \ 

UOOu    -    U)rs   y  1+  \  O-uJa  ) 

is  clearly  a  function  of  the  source-radius  and  the  cut-off  frequency 

for  the  mode  of  interest.   For  a  waveguide  with  pressure-release 

w  a 
boundaries  which  is  excited  in  the  lowest  mode,  the  quantity  

has  the  value  ~r  which  demonstrates  the  importance  of  the  ratio 

of  source-radius  to  layer-depth.  As  a  result  of  the  approximation 

which  led  to  Eq.  2.6,  Eq .  2.5  can  be  broken  into  two  parts: 


A.. 


^  [£&?!•**■> 
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Each  portion  contributes  significantly  to  u  only  when  the  exponent 

r 

of  each  is  stationary.   For  large  s  (which  corresponds  to  the 

leading  portion  of  the  signal),  only  the  first  integral  is  important 

and  yields  its  major  contribution  when  the  exponent  has  a  vanishing 

derivative  ""p 

1 


cAa- 


(a-T(j)-.O^  -     i     ======5    • 


Thus,  the  requirement  that  |  ho/ i  \    ~  U3as.se ts  an  upper  limit 


t*  --  T  liTtW 


on  the  time  within  which  the  first  integral  yields  the  major  contribution, 

For  a  step  input,  for  which  V(s)  =  — ,  the  radial  component  of  the 
particle  displacement  becomes 


Application  of  the  asymptotic  approximations  of  Eq.  2.6  leads  to 


Inversion    leads  to 


00  _ .  £*v\-* 


+A..fe(^^^r^^   27 

P2 3 

where  arg  =  w       J    t      -   T    . 

For  a  step  input,  the  solution  for  the  acoustic  pressure  becomes 
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The  quantity  t'  =  t  -  T  appearing  in  Eqs.  2.7  and  2.8  represents 
the  time  which  has  elapsed  since  the  propagating  transient  signal 
traveling  with  velocity  c  first  arrived  at  the  position  r.   This 
quantity  will  be  denoted  as  the  "delayed  time". 

For  the  leading  portion  of  the  propagating  transient,  one 
can  apply  a  standard  trigonometric  asymptotic  approximation 

to  the  functions  in  Eqs.  2.7  and  2.8.   Application  of  this  approx- 
imation to  the  sums  of  Eqs.  2.7  and  2.8  leads  to  expressions  of 

the   types  * 

__-    «-      imH         /       IT       -      Iff  \  f^    «  rJ        E     firH 


Application  of  Eqs.  2.9  and  2.10  to  the  solution  for  particle 
displacement  casts  it  in  the  form 

«»fa»-if-IJ,A»(£y,„an»-5-I0^    '■" 

aw 

where  2TTN  =  y  and  N  =  %w  T  — - 
J  a     n   c 

The  quantity  N  is  the  number  of  cycles  of  the  transient  generated 
in  time  t,  and  N  corresponds  to  the  cycle  where  the  approximations 

Cm 

(and  thus  the  theory)  fail  to  apply. 

By  combining  the  sine  and  cosine  terms  of  Eq.  2.11,  we  find 
the  expression  for  the  radial  component  of  particle  displacement 
to  be 
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2.12 


where  C(t)    -  ^  [  1 "  fo\(&.Y] *  t  [An^]1  and 

tan  ^    -  N°~  —     / 

1-  BnfNAlftV 

Application   of  Eqs .    2,9  and   2.10  to   the   solution   for   acoustic 
pressure   casts    it    in   the    form 


?(l)       K      W     IN      "^^^"       1     '   v   ^  2.13 


where  A(t)    =    j  Ql- Bo»  fa/N^]  *  +  C  A6i(N/^1  and 

tan  0  -   *"  — T .  .  , . ,  Ci    • 

t,   •       ,(5,14) 

Previous  work       has  shown  that  for  the  same  input  as  above 

the  exact  solution  for  propagation  of  a  transient  in  a  duct  becomes 
the  leading  term  in  the  solution  developed  for  the  isovelocity  layer 
It  should  be  noted  that  these,  duct  terms  can  be  retrieved  from  the 
trigonometric  approximations  by  letting  the  amplitude  factors  go 
to  unity  and  the  phase  angles  go  to  zero. 

The  solutions,  using  Bessel  functions  or  asymptotic  approxi- 
mations to  the  Bessel  functions,  for  the  transient  response  of  an 
isovelocity  layer  have  been  developed  with  the  restriction 

T  *.  t>  ±*~     or     M/NU  ±   i. 

Computer  programs  have  been  written  to  evaluate  the  predicted 
waveforms  for  acoustic  pressure  and  particle  displacement,  using 
either  the  Bessel  function  or  trigonometric  approximations  (see 
Appendices  A  and  B,  respectively).   These  programs  give  printed 
and  graphical  outputs  for  the  duct  and  layer  terms.   In  addition, 
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these  programs  predict  the  axis-crossings  and  amplitude-modulation 
in  the  layer.   Figures  2.1  and  2.2  illustrate  the  acoustic  pressure 
and  particle  displacement,  respectively,  predicted  by  the  two  approxi- 
mations.  Tables  2.1  and  2.2  give  the  calculated  values  of  the  predicted 
axis-crossings  and  amplitude-modulation  for  acoustic  pressure  and 
particle  displacement.   It  should  be  noted  that  the  first  portion 
of  the  first  half  -  cycle  of  the  waveform  predicted  by  the  trig 
computer  program  is  not  accurate.   To  avoid  a  mathematical  ambiguity 
at  the  leading  portion  of  the  waveform,  its  value  has  been  set  equal 
to  zero.   For  comparison  purposes,  the  predicted  amplitude  peaks 
were  normalized  to  the  value  of  the  highest  peak  in  the  waveform. 
As  can  be  seen  from  a  study  of  Tables  2.1  and  2.2,  the  difference 
between  the  two  predictions  is  on  the  order  of  one  per  cent  which 
is  smaller  than  the  experimental  accuracy  encountered.   Therefore, 
the  predictions  from  the  computer  program  using  trigonometric 
approximations  will  be  used  for  all  comparisons  with  experimental 
data. 

The  approximate  solutions  for  the  propagation  of  a  transient 
in  an  isovelocity  layer  have  been  developed,  and  computer  programs 
to  evaluate  these  solutions  and  predict  the  axis-crossings  and 
amplitude-modulation  in  the  layer  have  been  written. 

Transfer  Characteristics  of  Mylar  Transducers 

To  have  some  measure  of  confidence  in  the  received  waveform, 
it  was  decided  to  investigate  the  transfer  characteristics,  under 
steady- state  conditions,  for  Mylar  transducers.   This  investigation 
will  cover  the  radiation  impedance  at  the  surface  of  the  transmitter, 
r  =  a,  and  the  radiation  impedance  of  a  receiver. 
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Computed  using  Bessel  functions 
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Computed  using  trigonometric  approximations  to 
Bessel  functions 


Acoustic  Pressure  p  ts.  Delayed  Time  t1  as  Predicted  by  Computer 
Programs  for  the  Conditions  r=91.U  cm,  a=U.OO  cm  and  d-8.00  cm. 

Figure  2.1 
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Table  2, 

,1 

AXIS  CROSSINGS 

AXIS  CROSSINGS 

AMPLITUDES  OF 

AMPLITUDE 

USING  BES 

USING  TRIG 

MAXIMA,  BES 

MAXIMA,  Tl 

TIME  IN  USEC 

TIME  IN  USEC 

NORMALIZED 

normalize: 

1.000 

1.000 

1.45 

1.39 

0.632 

0.647 

7.53 

7.46 

0.470 

0.478 

18.32 

18.26 

0.390 

0.396 

33.58 

33.51 

0.339 

0.344 

52.99 

52.91 

0.303 

0.308 

76.17 

76.08 

0.275 

0.280 

102.72 

102.64 

0.253 

0.258 

132.33 

132.21 

0.235 

0.239 

164.57 

164.42 

0.219 

0.224 

199.13 

198.96 

0.205 

0.210 

235.72 

235.51 

0.193 

0.198 

274.07 

273.82 

0.182 

0.186 

313.96 

313.65 

0.171 

0.177 

355.18 

354.81 

Comparison  of  computed  axis-crossings  and  amplitudes  of  maxima  for 
acoustic  pressure  using  Bessel  functions  and  trigonometric  approxi- 
mations to  Bessel  functions  for  the  conditions  r  =  91.4  cm,  a  =  4.00  cm, 
and  d  -  8.00  cm. 
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Computed  using  trigonometric  approximations  to 
Bessel  functions 


Radial  Component  of  the  Particle  Displacement  E^  vs.  Delayed  Time  t1 
Predicted  by  Computer  Programs  for  the  Conditions  r«91.8  cm,  a=6.00  cm, 
and  d-8.00  cm. 

Figure  2.2 
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AXIS  CROSSINGS   AXIS  CROSSINGS 

USING  BES       USING  TRIG 
TIME  IN  USEC     TIME  IN  USEC 


3.75 

12.48 

25.94 

43.82 

65.75 

91.33 

120.16 

151.85 

186.05 

222.44 

260.73 

300.68 

342.07 


3.92 

12.60 

25.99 

43.79 

65.64 

91.13 

119.89 

151.51 

185.65 

221.99 

260.24 

300.15 

341.50 


Table 

2.2 

ilNGS 

AMPLITUDES   OF 

AMPLITUDES   OF 

G 

MAXIMA,    BES 

MAXIMA,    TRIG 

;ec 

NORMALIZED 

NORMALIZED 

0.474 

0.475 

0.701 

0.698 

0.844 

0.841 

0.933 

0.931 

0.982 

0.981 

1.000 

1.000 

0.995 

0.996 

0.974 

0.975 

0.941 

0.943 

0.901 

0.903 

0.857 

0.860 

0.811 

0.814 

0.764 

0.768 

Comparison  of  computed  axis-crossings  and  amplitudes  of  maxima 
for  the  acoustic  particle  displacement  using  Bessel  functions  and 
trigonometric  approximations  to  Bessel  functions  for  the  conditions 
r  =  91.8  cm,  a  =  6.00  cm,  and  d  *  8.00  cm. 
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In  a  steady-state  condition,  the  acoustic  potential  for  a 
cylindrical  source  in  a  waveguide  can  be  expressed  as 

<t>r--Kl)Z(r)ziMt.  2.14 

The  time  independent  wave  equation  becomes 

a'0r        I    diss*-      ,,z   . 

The  separation  of  variables  suggested  in  Eq.  2.14  leads  to  the 

w 
value  of  the  separation  constant,  k  =  — ,  where  the  phase  velocity 

CP 
c   is  defined  as 
P 

C 
Cft  -   rr — 77— vT  * 


t     {TWf 


2.15 


Thus,  the  wave  equation  for  the  radial  component  of  the  acoustic 
potential  becomes 

with  the  general  solution 

Rlr)--  A^^rf'liB.Yo^^.  2.i6 


For  an  outward  traveling  wave,  Eq.  2.16  becomes 

2.17 


RM'A^^^-'^^^V^ 


The  acoustic  potential,  acoustic  pressure,  and  the  radial  component 
of  the  particle  velocity  become 

*^     ,     u^y  u/r\  tot 
ST)  -   W>  (cr  )  (cf  )c   « 
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Using  Eq.  2.15,  the  specific  acoustic  impedance  Z_  becomes 

Hi  (<xva; 


2.18 


where  arg  =  «  J  ^  (J*ff    . 


Then,  the  radiation  impedance  Z   is  given  by 

Ir-ht  ~   i'pcfA    H°tt\^\  . 

where  A  is  the  total  radiation  area  of  the  transducer.   Thus,  the 
radiation  impedance  at  the  surface  of  the  transducer,  r  =  a,  is 
given  by 


The  behavior  of  Eq.  2.19  as  a  function  of  frequency  will  be 
investigated.   An  estimate  of  the  maximum  frequency  encountered 
in  a  shallow-water  layer  can  be  derived  by  recalling  that  the 
highest  frequencies  in  the  transient  occur  in  the  first  half-cycle. 
The  following  relation  concerning  the  time  for  the  maximum  frequency 
can  be  derived:  

mJlU-T1,    -ZT1M-7T    or 

The  relation  for  the  instantaneous  frequency 


U)  - 


can  then  be  manipulated  into  the  desired  result 

1  tC^)1,  • 
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For  r-a  =  100  cm  and  d  =  8  cm, 

w    =  650  khz. 
max 

The  cut-off  frequency  for  the  first  mode  and  the  above  approximations 


is 

o 
w  =59  khz. 
n 

Thus,  the  range  of  frequencies  expected  to  be  encountered  are 

w  £,  w  *  w    .   The  expected  frequency  must  always  be  greater  than 
n        niaX 

the  lower  limit,  for  the  assumption  of  a  steady  state  condition  no 
longer  holds  at  w  =  w  . 

Since  the  Hankel  functions  are  complex  numbers,  they  can  be 
expressed  as  a  modulus  and  a  phase  angle.  Thus,  Eq .  2.19  can  be 
rewritten  as 

where  Mod  is  the  modulus  and  9  is  the  phase  angle.   Figure  2.3 
shows  the  dependence  of  the  phase  angle  and  modulus  on  the  argument 
of  the  Hankel  functions.   For  small  values  of  the  arguments,  the 
radiation  impedance  is  strongly  dependent  on  the  frequency.   How- 
ever, when  the  value  of  the  argument  is  about  II  which  corresponds 
to  a  frequency  within  0.5  per  cent  of  the  cut-off  frequency,  the 
radiation  impedance  can  be  approximated  by  Ape  .   Thus,  for  the 
frequency  range  of  interest,  the  radiation  impedance  is  frequency 

independent , 

( 1  f\} 
The  open-circuit  mechanical  impedance  of  a  Mylar  transducer 

is  given  by  Z  =  r  +  i(wm  -  s/w)  where  r  is  the  transducer  mechanical 

°      J      m 

resistance  and  s  is  the  transducer  mechanical  stiffness.   For  the 
frequency  range  of  interest,  the  mechanical  resistance  is  considered 
to  be  negligible  in  comparison  with  the  reactive  part  of  Z  .   Under 
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(  1  fi^ 

the  assumption  that  the  mechanical  stiffness  is  given  by 

s  =  P  A/jj)  where  P  is  the  atmospheric  pressure  and  J(   is  the 

effective  depth  of  the  air  cushion  behind  the  Mylar  diaphragm,  the 

equation  for  Z  can  be  rewritten  as 

m  _. 

where   t   is   the   thickness   of   the  Mylar.     With  the  nominal  values 


r^mylar 


3 
=   1.395    gm/cm    , 


-3 

t   =  2.54  x   10        cm, 


P      =    1.013   x    106   dyne/cm    ,    and 
o 

-3 


J(=   2,00  x    10        cm, 


the  values  of  Z   for  the  frequency  range  of  interest  can  be 

calculated.   Figure  2.4  shows  the  dependence  of  the  magnitude  of 

Z 

m      . 
—  vs .  frequency . 

The  velocity  of  the  diaphragm  of  a  Mylar  transducer     is 

given  by 

[X*       C   V  ■£- ,,.,      2.20 

where  C   is  the  blocked  diaphragm  electrical  capacitance,  v  is 

the  signal  voltage  across  the  transducer,  0  is  the  turns  ratio, 

and  R  is  the  internal  resistance  of  the  power  amplifier.   The 

turns  ratio  0  is  defined  as 
_   CoVo 


0-- 


where  V   is  the  external  polarizing  voltage  and  x   is  the  position 
of  the  Mylar  film  when  the  polarizing  voltage  is  applied.   A 
comparison  of  the  magnitudes  of  the  terms  in  the  denominator  of 
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is  proportional  to  the  particle  displacement.   These  conclusions 
only  hold  for  the  leading  portions  of  the  received  transient  where 
Z   can  be  replaced  by  APc, 

The  effect  of  the  transfer  characteristics  of  Mylar  transducers 
has  been  found  to  be  dependent  on  the  frequency  and  the  radius  of 
the  transducer.   Given  a  transient  signal  which  is  the  mirror  of 
the  input  voltage,  the  following  is  a  summary  of  the  effects  thf.t 
would  be  observed  as  a  result  of  the  radiation  impedance: 

(a)  since  the  radii  of  the  family  of  transmitters  were  large, 
the  response  is  essentially  flat  throughout  the  region  of  interest; 
thus,  no  effects  should  be  observed. 

(b)  because  of  the  small  radius  of  the  receiver,  the  effect 
becomes  more  pronounced  than  for  the  transmitters.   One  should 
observe  a  reduction  in  the  size  of  the  peaks  of  the  leading  portion 
of  the  transient  and  should  observe  an  enhancement  of  the  peaks  as 
the  local  apparent  frequency  approaches  the  cut-off  frequency.   A 
similar  phase  angle  effect  should  be  observed. 
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3.   EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  in  this  research  project  consisted 
of  a  tank,  a  family  of  cylindrical  Mylar  transducers,  a  barium 
titanate  array  receiver,  and  the  associated  electronics  equipment. 
The  laboratory  setup  is  shown  in  Fig.  3.1. 

Tank 

The  water  layer  was  contained  in  a  tank  constructed  of  5/8 
inch  marine  plywood.   The  dimensions  of  the  tank  were  6  ft  x  6  ft  x 
1  ft.   The  plywood  was  covered  with  Varathene  to  provide  water- 
proofing and  lined  with  1/4  inch  polyethylene  foam  which  was  assumed 
to  provide  the  desired  pressure-release  bottom  boundary  of  the 
layer.   The  polyethylene,  not  being  waterproof,  was  coated  with 
liquid  Neoprene  paint  to  prevent  deterioration  of  the  pressure- 
release  surface.   The  bottom  of  the  water  layer  was  maintained 
level  by  means  of  adjustable  feet  installed  in  the  legs  of  the 
supporting  platform.   Since  the  walls  of  the  tank  were  also  lined 
with  polyethylene,  the  tank  could  be  excited  as  a  pressure-release 
cavity,  which  provided  one  method  of  determining  the  depth  of  the 
layer. 

Transducers 

Previous  work   '    had  indicated  that  Mylar  dielectric 
transducers  were  well-suited  for  use  in  layers. 

A  family  of  eight  Mylar  transducers  and  one  barium  titanate 
linear  array  were  constructed  for  use  in  this  investigation.   The 
radii  of  the  Mylar  transducers  were  10  cm,  6  cm,  4  cm,  2  cm,  1  cm, 
0.5  cm,  and  two  at  0.25  cm.   Each  had  an  active  radiating  face 
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8  cm  in  height  which  was  the  maximum  depth  of  the  layer.   These 

various  transducer  dimensions  were  chosen  to  provide  a  large  range 

of  values  for  the  ratio  N/N  .   The.  barium  titanate  array  was  8  cm 

a 

high  and  1/8  inch  in  diameter. 

The  Mylar  transducers  were  constructed  using  techniques 
described  in     ,  modified  as  follows: 

(a)  The  backplates  of  the  larger  transducers  were  constructed 
from  \   inch  thick  aluminum  cylinders.   The  backplates  for  the 
smallest  transducers  were  fabricated  from  \   inch  diameter  copper 
tubing. 

(b)  The  cylinders  were,  capped  at  the  top  and  bottom  by 
lucite  plates. 

(c)  The  Mylar  film  was  cut  to  the  desired  dimensions  and 
attached  to  the  aluminum  cylinders  with  thin  strips  of  double-stick 
tape,  one  piece  around  the  top,  one  around  the  bottom,  and  one 
piece  vertically  from  bottom  to  top. 

(d)  The  inner  electrical  connection  was  made  from  the  center 
pin  of  a  BNC  connector  mounted  in  the  center  of  the  top  lucite 
plate  to  the  inside  of  the  backplate.   The  outer  conducting  path 
from  the  metallized  side  of  the  Mylar  film  to  the  BNC  was  completed 
by  covering  the  top  lucite  plug  with  aluminum  foil  which  was  lapped 
over  the  Mylar  film.   To  ensure  good  electrical  continuity,  silver 
paint  was  applied  to  the  aluminum-Mylar  interface.   To  decrease  the 
possibility  of  electrical  pickup,  the  bottom  lucite  plate  was  also 
covered  with  aluminum  foil. 

(e)  When  the  silver  paint  had  dried,  the  entire  transducer 
was  covered  with  liquid  neoprene  to  provide  the  necessary  water- 
proofing. 
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Tank  with  transmitter  and  receiver.  Electronics  equipment  in 
the  background, 

laboratory  Setup 

Figure  3.1 


Transmitters  and  Receivers. 
re  3.2 
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The  pressure-sensitive  barium  titanate  was  constructed  using 
eight  1/8  inch  by  1/8  inch  cylinders  mounted  on  several  strands 
of  twisted  wire,  which  provided  the  inner  electrical  connection. 
To  provide  the  damping  necessary  to  reduce  the  effect  of  longitudinal 
resonances,  the  spaces  between  the  cylinders  were  filled  with  strips 
of  rubber.   The  outer  electrical  connection  was  provided  by  a  wire 
which  was  coiled  around  each  cylinder.   This  probe  proved  to  be 
somewhat  insensitive  and  was  used  only  for  comparison  purposes. 

When  the  electrical  load  on  the  Mylar  receiver  (see  Fig.  2.5) 
was  terminated  by  the  1  megohm  impedance  of  the  voltage  amplifier, 
the  receiver  was  sensitive  to  particle  displacement.   Therefore, 
the  Mylar  receivers  were  used  to  measure  the  particle  displacement, 
and  the  barium  titanate  array  was  used  to  measure  the  acoustic 
pressure.   Transmitters  and  receivers  are  shown  in  Fig.  3.2. 

Electronic  Equipment 

Figure  3.3  shows  the  block  diagram  of  the  equipment;  Table 
3.1  identifies  the  components  in  the  figure.   The  unit-pulse  generator 
was  used  for  the  step  input. 

The  received  signal  was  sent  from  the  receiver  through  the  load 
and  bias  network,  amplifier  and  filters  to  trace  "B"  of  the  dual- 
beam  oscilloscope.   The  output  of  the  power  amplifier  could  be 
displayed  on  trace  "A".   Measurements  of  the  received  waveform  axis 
crossings  and  amplitudes  were  made  directly  from  the  oscilloscope 
display,  or  photographed  for  later  analysis.   The  electronics  were 
triggered  by  the  time-mark  generator  which  also  triggered  the  time- 
base  of  the  oscilloscope.   The  intensity  of  the  "B"  trace  was 
modulated  by  a  10  usee  marker  available  from  the  time-mark  generator. 
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Table  3.1 
Electronics  Equipment 

UPG  Unit  pulse  generator,  General  Racio  type  1217-B 

PWR  AMP  Power  amplifier,  Hewlett  Packard  467 -A 

TMG  Time  Mark  Generator,  Tektronic  type  180A 

AMP  1  Voltage   amplifier,  Tektronic  type  112 1A 

SKL  FIL  Electronic  band  pass  filter,  SKL  model  202-D 

RC  FILl  Homemade  variable  low-pass  RC  filter  with  400  khz  roll-over 

RC  FIL2  Homemade  high-pass  RC  filter  with  2  khz  roll-over 

SCOPE  Dual  beam  oscilloscope,  Fairchild  model  777 

LBN  Load  and  bias  network  shown  in  Fig.  3.4. 
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Source  and  receiver  were  carefully  positioned  so  that  the  distance 
between  adjacent  sides,  r-a,  was  always  the  same..   This  insured 
that  the  trace  modulation  would  be  properly  positioned  with 
respect  to  the  inception  of  the  received  transient  signal.   The 
load  and  bias  network  shown  in  Fig.  3.4  served  to  supply  the  D.C. 
voltage  bias  needed  by  the  Mylar  transducers  without  allowing 
this  voltage  to  appear  elsewhere  in  the  circuit  and  to  provide 
the  proper  electrical  load  for  the  receiver.   The  22  ohm  resistor 
in  the  transmitter  circuit  was  needed  to  suppress  parasitic 
oscillations  caused  by  the  transient  excitation  of  the  capacitance 
of  the  transmitter  and  the  residual  inductance  of  the  power 
amplifier. 

The  rise  time  of  the  system  was  experimentally  determined 
to  be  about  0.5  usee. 
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4.   EXPERIMENTAL  INVESTIGATION 

The  purposes  of  the  experimental  investigation  were  to   (a)  observe 
and  analyze  the  received  acoustic  pressure  and  radial  particle  displace- 
ment waveforms  in  the  layer,   (b)  compare  these  waveforms  with  the 
theoretical  predicted  waveforms  produced  by  the  computer  programs,  and 
(c)  vary  the  value  of  the  ratio  N/N  over  a  large  range  to  investigate 

9, 

the  validity  of  the  assumed  asymptotic  approximations  previously  dis- 
cussed.  Four  basic  parameters  characterized  each  experiment:   the 
range  r,  the  transmitter  radius  a,  the  layer  depth  d,  and  the  speed 
of  sound,  c. 

Experimental  Procedures 
Before  each  run, it  was  decided  which  parameter  would  be  varied. 
Generally,  the  time  of  flight  T  and  the  depth  of  the  layer  d  were 
kept  constant  while  the  source  radius  a  was  varied  throughout  its 
range.   The  0.25  cm  radius  Mylar  transducer  was  used  to  record  the 
received  particle  displacement,  and  the  barium  titanate  linear  array 
was  used  for  the  acoustic  pressure  receiver.   The  temperature  of  the 
water  was  measured  to  determine  the  speed  of  sound  in  the  layer 
The  depth  was  determined  by  averaging  several  physical  measurements 
taken  over  the  ensonified  area.  Another  method  of  determining  the 
depth  was  to  excite  the  tank  as  a  pressure-release  cavity.   The 
cut-off  frequency  could  then  be  determined  by  observing  the  received 
signal  on  an  oscilloscope  and  varying  the  frequency  until  the  signal 
disappeared.   Once  the  cut-off  frequency  was  known,  the  depth  could 


TT 


c 


be  calculated,  using  the  formula  w  =  — r  .   Because  of  the  irregularities 


n 


observed  in  the  depth,  this  method  proved  to  be  no  more  accurate  than  the 
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actual  measurement  using  a  ruler  and  was  consequently  abandoned. 
The  intensity  of  the  trace  of  the  oscilloscope  was  modulated  by 
a  10  usee  marker  available  from  the  time-mark  generator.   The 
source  and  receiver  w^re  carefully  positioned  so  that  the  distance 
between  adjacent  sides,  r-a,  was  always  the  same  and  such  that 
the  reception  of  the  received  transient  and  one  of  the  time  marks 
always  coincided.  This  insured  that  when  changing  transmitters, 
identical  positioning  and  time  of  flight  could  be  achieved. 

The  received  waveforms  were  photographed  from  the  oscil- 
loscope display.   For  the  stap  input,  several  photographs  were 
taken  of  each  waveform.   An  expanded  time  scale  was  used  when 
photographing  the  leading  portion  of  the  waveform  which  contained 
the  high  frequencies.   For  a  typical  experimental  run,  three 
photographs  of  each  waveform,  were  made  utilizing,  for  -xample, 
20  usec/cm  time  scale.   By  varying  the  delay  time  vernier,  the 
entire  waveform  could  be  displayed  an i   photographed  for  later 
analysis.   By  careful  measurement  of  the  position  of  e      sis 
crossing  with  the  neighboring  time  marks,  the  relative  phase  of 
each  axis  crossing  could  be  compared  between  the  photographs  taken 
for  different  source  radii.   The  amplitudes  of  the  peaks  of  the 
waveforms  were  read  with  the  help  of  dividers  and  a  seal-.   Sample 
photographs  are  shown  in  Fig.  4.1.   It  should  be  notid  that  the 
anomalous  behavior  of  the  received  waveforms  as  evidenced  by 
Fig.  4.1  departs  drastically  from  the  smoothly-varying  computer 
predictions  of  Figs.  2.1  and  2.2.   This  anomaly  will  be  discussed 
in  more  detail  later. 
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For  each  experimental  run,  a  corresponding  computer  run  was 
made  using  the  same  values  for  r,  a,  c,  and  d.   The  computer 
programs  wer^.  written  so  that  their  output  f      h  run  includ. 
graphs  of  the  predicted  waveforms,  times  of  axis  crossings  for 

the  layer  solution  t1  and  for  the  duct  solution  t1  ,  and  amplitudes 

m  o" 

of  the  peaks  and  the  time  of  occurrence  of  these  peaks  for  the 
layer  solution.   In  this  manner,  the  theoretical  values  of  the 
times  of  axis  crossings  and  of  the  amplitude  peaks  could  be 
compared  with  thos==  values  experimentally  determined. 

The  time  of  the  axis  crossings  »ere  compared  in  the  following 
manner:   Previous  work    suggested  that  the  time  of  aach  axis 
crossing  should  be  the  same  as  the  theoretically  predicted  t' 
within  the  validity  of  the  apprc       ins  us  d.     «  s  ^hown  in 
section  2,  the  theory  was  valid  only  as  long  as  N/N  £  1.   There- 

■x 

fore,  the  convenient  parameters  for  analyses  were  the  time  of  the 

axis  crossings  and  the  ratio  N/N  .   In  ordt-.r  to  reduce  the  size  of 

a 

the  graphs,  the  time  coordinate  was  reduced  by  th^  amount  t'  predicted 
by  the  computer..  The  time  coordinate  signified  by£\t'  was  plotted 
against  the  ratio  N/N  .   Figuces  4.2  to  4.4  are  representative 

EL 

examples  of  these  graphs.   Thp  solid  line  represents  the  values  of 

t1   -  t'f  and  the  circled  dots  represent  the  values  of  t1    -  t1. 
om    m  exp    m 

Were  the  theory  to  have  been  completely  correct  and  the  experiment 

accurate  enough,  the  data  would  have  fallen  on  the  X-axis  (which 

represents  t1  -  t').   As  can  be  seen  from  Figs.  4.2  to  4.4,  the 
mm 

data  were  not  in  agreement  with  the  theory.   The  effect  of  the 
phase  angle  demonstrated  in  Fig.  2.6  was  taken  into  account  and 
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plotted  on  Figs.  4.3  and  4.4.   As  can  be  seen,  this  did  not 
improve  the  agreement  between  the  theory  and  the  data. 

The  peak-amplitude  data  were  also  analyzed.   The  experimental 
amplitude  peaks  were  first  normalized  to  tht  value  of  the  highest 
peak  in  the  waveform.   The  predicted  peaks  were  normalized  in  tl 
course  of  the  computer  programs.   The  ratio  of  these  normalized 
quantities  was  then  plotted  as  a  function  of  the  delayed  time  of 
the  predicted  amplitude  peaks  for  each  run.   Exact  agreement  in 
the  shape  of  the  transient  would  result  in  a  lin~  of  zaro  slopt 
passing  through  the  point  (1,0).   The  amplitude  effect  discussed 
in  Section  2  for  a  Mylar  transducer  used  as  a  receiver  was  also 
plotted.   Figures  4.5  to  4.7  show  plots  for  various  values  of  r, 
a,  and  d.   Once  again,  this  affect  did  not  improve  the  agreement 
between  the  theory  and  the  data. 

Error  Analysis 

The  rise  time  of  the.  experimental  system  was  determined  to  be 
about  0.5  usee.   Therefore,  instead  of  b^-ing  a  perfect  step- function 
input,  the  velocity  at  the  source  was  initially  described  by 

r  -I 

where  o\=  2  x  10   sec  '  is  the  reciprocal  of  the  rise  time. 
Previous  work    has  indicated  that  the  rise  time  of  the  system 
should  introduce  distortion  into  the  received  signal  for  at  most 
3.0  usee  of  delayed  time  t'.   Since  this  time  is  less  than  that 
for  the  first  maxima,  the  effect  of  the  rise  time  on  the  received 
signal  can  be  ignored  in  this  investigation. 
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It  has  been  shown  that  the  velocity  of  the  Mylar  face  decays 
to  88  per  cent  of  its  maximum  value  after  a  time  of  360  usee 
This  effect  is  not  large  enough  to  e.ccount  for  the  tailing  off 
of  the  amplitude  peaks  of  the  received  waveforms  sern  in  Figs. 
4.5  to  4.7. 

Uncertainties  in  r,  a,  c,  and  d  were  the  possible  sources 
of  error  in  this  investigation.   The  effecti    r  -^ius  of  each 
transmitter  was  measured  to  the  nearest  0.5  mm,  which  could  have 
introduced  an  uncertainty  of  +  1.0  per  c-  nt  in  the  value  of  a 
for  the  6.0  cm  radius  source  and  +5.0  per  cent  in  the  value  of  a 
for  the  1.0  cm  radius  source.   Measurement  of  the  temperature  of 
the  water  to  the  nearest  half-degree,  restricted  the  systematic 
error  in  c  to  within  +  0.5  per  cent.   Determi nation  of  th^  tirr.e- 
of- flight  T  of  the  received  signal  form  the  source  to  the  receiver 
produced  a  possible  systematic  error  of  +  0.5  p.-r  cent,  which  in 
turn  resulted  in  an  error  in  r  of  +  0.5  per  cent. 

The  largest  experimental  error  was  encountered  in  determining 

the  depth  of  the  layer.   Observed  irregularities  of  +  1.0  mm  in 

the  bottom  surface  of  the  tank  could  introduce  an  error  in  the 

local  effective  depth  of  the  layer.   Since  the  averaged  layer  depth 

used  as  an  input  to  the  computer  programs  was  only  accurate  to  +  1.0 

mm,  the  systematic  error  introduced  in  w  was  on  the  order  of  +  3.0 
»      j  n  _ 

per  cent  for  a  depth  of  8.0  cm. 

The  error  in  determining  the  times  of  the  axis  crossings  from 
the  data  was  effected  by  the  accuracy  of  the  oscilloscope  time-base 
which  was  calibrated  to  within  0.3  per  cent.   With  the  trace  of 
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the  oscilloscope  modulated  with  the  10  usee  marker  as  an  additional 
time  reference,  this  error  was  reduced  even  more.   Therefore,  a 
random  error  of  +  1.0  per  cent  could  be  expected  when  one  takes 
into  account  the  accuracy  of  reading  the  oscilloscope  scale. 
The  errors  in  the  computed  values  of  t'  were  influenced 


IT! 


predominately  by  the  error  in  w  ,  since  it  appears  in  the  argument 

(2    2s 
of  the  Bessel  functions  (w  J  t   -  T   )  from  which  the  t'  are 

n  V  m 

(5) 
computed.   It  has  been  shown    that  the  argument  of  the  Bessel 

functions  could  be  rewritten  as 


"•ntUtWD    =•  C^ 


where  C  are  a  set  of  constants  which  yield  zero  values  for  the 
m 

amplitude  of  the  particle-displacement  solution.   For  the  leading 
portion  of  the  signal,  we.  have 

w*(zT  Ct'mV2-  «  Cm 

and  C,^ 


Therefore,  the  +3.0  per  cent  error  in  w   could  lead  to  a  +  6.0 
—  n  — 

per  cent  error  in  the  computed  values  of  t'. 

m 

To  reduce  the  scatter  of  the  data,  a  one-parameter  fit  was 
tried.   It  was  decided  to  determine  an  effective  cut-off  frequency 
and  use  this  as  an  input  to  the  computer  programs  instead  of  the 
layer  depth.   This  was  accomplished  by  determining  the  cut-off 
frequency  for  the  values  of  t^   and  the  arguments  of  the  trig 
function  approximations  developed  in  Section  2.  At  an  axis 
crossing,  the  argument  was  equal  to  an  integer  multiple  of 


[ 
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A  computer  program  was  written  to  evaluate  the  cut-off  frequency 
for  each  axis  crossing.   Then,  these  cut-off  frequencies  were 
averaged  to  give  an  "effective  cut-off  frequency"  to  be  used  in 
the  computer  programs  to  recompute  the  time  for  each  axis  crossing. 
Figure  4.8  shows  the  effect  of  this  one-parameter  fit.   As  can  be 
seen,  it  did  not  improve  the  scattering  significantly  and  was  not 
used. 

Error  flags  on  Figs.  4.2  to  4.8  indicate  the  relative  magnitudes 
of  errors  in  the  experimental  and  computed  values  of  the  axis  crossings 
and  peak  amplitudes.   However,  these  error  flags  are  not  large  enough 
to  explain  the  scatter  of  the  data  or  the  anomaly  which  was  present 
(see  Fig.  4.1).   The  theory  could  not  account  for  this  anomaly  so 
various  experiments  were  conducted  to  investigate  its  character  and 
source  of  origin. 

Firstly,  it  was  suggested  that  irregularities  in  the  bottom  of 
the  tank  might  be  a  possible  cause  of  the  anomaly.   Therefore,  the 
relative  positions  of  the  transmitter  and  the  receiver  in  the  tank 
were  changed.   The  values  of  T  and  d  were  always  kept  constant.   Four 
contrasting  positions  were  chosen;  one  near  the  left  hand  side  of 
the  tank,  one  in  the  center  of  the  tank,  one  adjacent  to  a  seam  in 
the  bottom  near  the  right  hand  side,  and  one  where  the  transmitter 
and  receiver  straddled  the  seam.   As  can  be  seen  in  Fig,  4.9,  within 
the  experimental  accuracy  of  the  investigation,  positioning  on  the 
bottom  had  no  effect  on  the  received  waveform. 

Secondly,  it  was  decided  to  monitor  the  current  drawn  by  the 
transmitter  when  loaded  by  the  water  medium  and  compare  this  with 
the  current  predicted  assuming  pc  loading.   The  current  was 
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experimentally  determined  by  measuring  the  current  in  the  two  legs 
of  the  circuit  shown  in  Fig.  4.10.   The  current  present  in  each 
leg  was  amplified  by  a  differential  amplifier  and  then  the  difference 
displayed  on  the  oscilloscope.   Photographs  of  the  received  signal 
could  then  be  analyzed  to  determine  the  current  drawn  as  a  result 

of  acoustic  radiation  into  the  medium.   By  analyzing  the  equivalent 

( 1  f>} 
circuit  of  a  Mylar  transducer  used  as  a  transmitter     ,  as  shown 

in  Fig.  4.11,  the  predicted  value  of  the  current  drawn  as  the  result 
of  the  assumed  Dc  loading  can  be  determined.   For  this  assumed 
loading,  the  current  drawn  by  the  mechanical  leg  of  the  source 
should  be  flat.   Typical  values  for  the  6  cm  radius  source  were 
near  6  uamp.   The  measured  values  of  the  current  had  rapid  oscil- 
lations of  nearly  1  uamp  peak-to-peak  amplitude,  which  indicated  a 
significant  fluctuation  from  the  flat  response  predicted.   This 
casts  doubts  on  the  validity  of  the  assumption  of  Oc  loading  used 
in  the  developments  of  Section  2. 

Lastly,  it  was  suggested  that  the  polyethylene  foam  was  not  a 
perfect  pressure-release  surface  and  that  transmission  through  it 
into  the  wood  bottom  was  the  cause  of  the  anomaly.   In  order  to 
investigate  this  possibility,  two  additional  layers  of  polyethylene 
were  glued  to  the  bottom  of  half  the  tank.   Then  runs  were  made  on 
the  single- layered  bottom  and  then  on  the  three- layered  bottom. 
The  quantities   r,  a,  and  d  were  kept  constant  for  the  corresponding 
runs.   Figures  4.12  to  4.14  indicate  the  predicted  waveform,  the 
experimentally  determined  waveform  on  the  three- layered  bottom,  and 
the  experimentally  determined  waveform  on  the  single- layered  bottom 
for  three  values  of  a:   10,  6,  and  4  cm.   The  vertical  and  horizontal 
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(a)  Predicted  Waveform 


Step  Input 
in  Velocity 

r=70.66  cm 
a=10.00  cm 
d=8.00  cm 


(b)  Received  particle  displacement  waveform  on  three-layered 
bottom 

Step  Input 
in  Velocity 

r=70.66  cm 
a=10.00  cm 
d=8.00  cm 


(c)  Received  particle  displacement  waveform  on  single-layered 
bottom 

Bottom  Layer  Investigation  for  Three-Layered  Bottom  and  Single-Layered 
Bottom  using  10,00  cm  Source. 

Figure  U.12 
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(a)  Predicted  Waveform 
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(b)  Received  particle  displacement  waveform  on  three-layered 
bottom 
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d=8.00  cm 


(c)  Received  particle  displacement  waveform  on  single-layered 
bottom 

Bottom  Layer  Investigation  for  Three-Layered  Bottom  and  Single-Layered 
Bottom  using  6,00  cm  Source. 

Figure  U.13 
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in  Velocity 


r=6U.66  cm 
a=li.00  cm 
d=6.00  cm 


(c)  Received  particle  displacement  waveform  on  single-layered 
bottom 

Bottom  Layer  Investigation  for  Three-Layered  Bottom  and  Single-Layered 
Bottom  using  U.OO  cm  Source. 

Figure  U.lU 


65 


scale  settings  on  the  oscilloscope  were  kept  the  same  for  each 
set  of  runs.   As  can  be  seen,  there  is  some  evidence  of  a  bottom 
effect,  but  there  is  no  clear  correlation  between  the  three  runs. 
However, this  evidence  is  sufficient  to  doubt  the  credibility  of 
the  assumption  of  the  perfect  pressure-release  characteristics  of 
the  0.25  inch  polyethylene  foam. 
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5.   RESULTS  AND  CONCLUSIONS 

The  investigation  has  shown  that  for  N/N  £  1  the  predicted 
and  experimentally  observed  shifts  in  axis  crossings  should  agree 
to  within  +  6  per  cent  if  one  takes  into  account  just  the  known 
systematic  errors.   The  correlation  between  the  predicted  and 
observed  amplitudes  should  be  of  the  same  accuracy.   However, 
the  difference  between  the  theory  and  data  is  greater  than  that 
expected. 

The  theory  as  developed  in  Section  2  could  not  account  for 
the  anomaly  which  was  present  in  all  runs,  but  the  qualitative 
agreement  between  the  observed  and  predicted  waveforms  and  the 
general  agreement  with  the  solutions  developed  in  other  work 
lead  one  to  believe  that  the  theory  presented  in  Section  2  is  correct. 

It  is  concluded  that  the  experimental  setup  used  is  not  accurate 
enough  to  give  results  needed  to  confirm  the  theory. 

It  is  recommended  that  better  electronic  equipment  and  transducers 
should  be  investigated.   Further  studies  should  be  conducted  to  check 
the  validity  of  the  assumed  pressure-release  characteristics  of  the 
polyethylene  foam,   More  investigations  are  needed  to  try  to  determine 
the  cause  of  the  observed  anomaly. 


(4,5,10,13) 
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